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Abstract 

A recently proposed modified perturbative expansion for QCD incorporating gluon condensation 
is employed to evaluate the quark and gluon self-energy corrections in the first approximations. 
The results predict mass values of 1/3 of the nucleon mass for the light quarks u, d, and s and a 
monotonously growing variation with the current mass. The only phenomenological input consists 
in that (G 2 ) is evaluated up to order g 2 as a function of the unique parameter C defining the 
modified propagator and then C is fixed to give a current estimate of (g 2 G 2 ). The light quarks u and 
d result to be confined and the s, c, b and t ones show damped propagation modes, then suggesting 
a model for the large differences in stability between the nucleons and the higher resonances. 
The above properties of quark modes diverge from the fully confinement result following from the 
similar gluon propagator previously considered by Munczek and Nemirovski. On another hand, 
the condensate effects on the gluon self-energy furnish a tachionic mass shell as predicted by the 
Fukuda analysis of gluon condensation in QCD. 

PACS numbers: 12.38.-t, 12.38.Bx, 12.40.Yx, 12.38.Aw, 12.15.Ff 
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I. INTRODUCTION 



One of the great achievements of theoretical High Energy Physics in the last thirty years 
has been the discovering and development of the Quantum Chromodynamics (QCD). The 
smallness of the coupling constant at high momenta (asymptotic freedom) has allowed to 
develop a perturbative framework applicable for very high energy processes. That situation 
strongly simplifies the study of such phenomena and the calculated quantities are in good 
agreement with the experimental data. However, perturbative QCD (PQCD) is far from 
being able to furnish even a rough description for the relevant physics at low energies. The 
solution of this situation is currently one of the main challenges of Particle Physics. 

A relevant phenomenon related with QCD is the color confinement. Nowadays, there are 
strong reasons to believe that the relation between the basic quantities in QCD, like the 
gluon and quarks fields, and the real world characterized by a whole variety of interacting 
mesons and baryons, can be understood by solving the confinement problem. The basic 
picture which seems plausible to be derived from QCD is that the fundamental fields can 
not be associated to asymptotic states of quark and gluons and then the true physical states 
consist only of colorless composites of quarks and gluons (mesons and baryons). 

The above quoted limitations of the PQCD mean in particular that the usual Fock space 
vacuum of the non-interacting theory, is unable to predict, even approximately, the real 
ground state properties of QCD |]1], ^|, ^ |]. This is at contrast with the case in QED 
where the standard perturbation theory gives a better than good concordance with the 
experimental data. Then, the nature of vacuum structure is one of the main problems to be 
clarified and naturally its solution is closely linked with the color confinement effect. Reviews 
of various models considered to investigate the confinement and the vacuum structure in 
QCD can be found in Refs. |5], ||. 

Another basic problem for QCD is the one related with the quark masses, that is, the 
so called current quarks masses which appear in the Lagrangian. As the quarks are not 
observed as free states, the meaning of the quark masses need to be considered with care 
|7], |8|, ^|. The nature of these parameters is an important object of research since long 
time ago |H| [11, 0. The relations among these masses and their magnitude have been 
determined through the methods of current algebra in combination with measurements of 
the light mesons masses. 
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In the former works fl3| , [14], |15| an attempt directed to construct a modified perturbation 
expansion for QCD being able to predict at least some low energy properties of this theory 
have been considered. First, in the modified expansion conserving the color SU(3) and 
Lorentz symmetries was proposed, aiming to solve the symmetry limitations of the earlier 
chromo-magnetic field models |TJ, |3|, |3J] which inspired the search. A similar lack of manifest 
Lorentz invariance had also a modification of the Feynman rules, intending to include gluon 



condensation, advanced later in flTJ| [17|, [L8[ in which a delta function for k < A qc< i was 
summed to the perturbative piece. The expansion proposed in considers a change of 
the gluon propagator in a term associated to a condensation of zero momentum gluons. This 
alteration of the usual rules had the interesting property of producing a non vanishing value 
for the gluon condensation parameter (g 2 G 2 ) already in the first approximations JT^]. In 



addition, in |13| a non vanishing value for the effective self-energy of the gluons was obtained 
at the loop level. Finally, a perturbative evaluation of the effective potential as a function of 
the condensate parameter indicated that the condensation is spontaneously generated from 
a zero condensate state. 



The next work [14, 15], had as a main objective to justify the applicability of the Feynman 
expansion introduced in ||13|| . This search was in need because after modifying the propaga- 
tor, it was unknown whether or not the initial wavefunction generating the expansion was 
a physical state of the free theory. In [TJ|, |T5| , by making use of the operational formulation 
of QCD developed by Kugo and Ojima PC|| , it was possible to find a physical state of the 
non-interacting theory being able to generate through the Wick expansion, the sort of prop- 
agators considered in The discussion allowed also a more precise characterization of 
the class of changes admitted in the diagrammatic expansion by the physical state condition 
on the initial state. Specifically, it was found that the C parameter, describing the gluon 
condensation, must be real and positive. 

The present work is an expanded version of a letter sent to publication. Its main objective 
is to evaluate first (G 2 ) in terms of the condensate parameter C up to order g 2 , in order 
to determine the quark masses after fixing the parameter g 2 C to reproduce the accepted 
value for (g 2 G 2 ). The essential result is the evaluation of the constituent masses for the light 
quarks. In order to perform it, we have made use of the recent values reported for the current 
masses The idea inspiring the work was the possibility that the gluon condensation 

described by modified expansion could predict the appearance of large masses for the light 
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quarks (u, d and s) of the order of 1/3 of the nucleoli mass (constituent quark mass), thus 
furnishing an explanation for the successful constituent quark models. Surprisingly, after 
fixing the condensate parameter g 2 C to reproduce the most accepted value of the (g 2 G 2 ) 
(through the here performed evaluation of this quantity in terms of g 2 C), it followed that the 
light quarks u, d and s got exactly the needed constituent like values for the masses in one of 
the propagation modes of the modified mass shell. The free propagation modes, which are at 
discrepancy with the absence of free quarks, are here only arising within an approximation 
evaluating the main effects of the condensate. After the inclusion of all the terms completing 
the one loop approximation in the gluon self-energy an interesting picture arises. It follows 
that the previously free propagation modes of the u and d quarks turn to be confined. On 
another hand, the s, c, b and t quarks waves become damped. That is, the squared rest mass 
get complex values. It should be stressed that the peaked structure of the propagator for the 
u and d at ^pp 2 = m = 0.30915 GeV signals the value of m as the one expected to appear 
in the kinetic energy terms "p* 2 / (2m) defining the free part of the bound state equations for 
u and d constituted nucleons. Henceforth, a possible mechanism contributing to justify the 
drastic difference between the stability of nucleons and the higher resonances is suggested. 
The absence of propagating modes for the u and d, might be assuring the absolute stability 
shown by the nucleons bound states. In other manner, the presence of damped modes in 
the rest of the quarks would contribute to the decaying channels observed for all the hadron 
resonances associated to these short lived states. 

Next, after considering the ground state within each of the baryon resonance groups as 



classified in |yj and its set of three constituent quarks, and adding the calculated quark 
masses for each of them (which depends on the flavor through the assumed values of the 
current masses), a spectrum of baryon masses was estimated. The results again reproduce 



reasonably well the experimentally determined resonance energies pl| . This outcome sug- 
gests the less relevant role of the interaction energy among quarks in determining the masses 
of these hadrons. 



At this point it is important to make clear that earlier works ^3] introduced a pure 
delta function at zero momentum in the gluon propagator, searching for a model of meson 
resonances. In these treatments the constants introduced, multiplying the delta function, 
were different to the one fixed in our previous work |TJ], (I5| , leading also to very different 
results. In particular, the singularity structure of the quark propagator have no poles on the 
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real p 2 axis J2ij at difference with our results. This difference is directly connected with the 
definite sign for the coefficient of the (^-function which in our case was precisely determined 
H| , [T5|| . Since our analysis in that paper was based in the direct construction of the 



m 



state incorporating the condensate, it follows that the difference with the results of E2, 
is rooted in a corresponding difference in the initial state determining the Wick expansion. 
It could be also possible that both procedures lead to the same exact physical predictions, 
but the first approximations in each approach could turn to be more appropriate for the 
discussion of specific physical issues. This question will be addressed elsewhere. It should 
be also further noticed that gluon condensate states reproducing the modified expansion of 



22, 23 have been also discussed in 24, 25 . However, all these references are related with 



states showing confined quarks and massive gluons, at variance with our proposed state 
determining massive quarks (in the simplest approximation) and tachyonic gluons. 

It is also following that another dispersion relation is obtained which predicts a vanishing 
value of the mass as tuq — > 0. This mode could be connected with the family of low mass 
mesons (e.g. II and K mesons) . Specifically, the possibility is open that the bound states of 
quark excitations in these light mass states of quarks could describe such low lying mesons 
within the here presented framework. A circumstance pointing in this direction is that the 
spontaneous chirality breaking associated with this mass generation should be expected to 
be associated with some massless Goldstone modes. Therefore, it is possible that the bound 
states of such modes predicted by the Bethe-Salpeter equations turn to describe the light 
mesons II and K. This question however, will be deferred to a further study. 

The calculations of (G 2 ) and the quark and gluon masses will be performed up order g 2 . 
It seems worth anticipating that the approximation scheme will resemble a sort of quasi- 
classical limit in which the condensate seems to play a role of a macroscopic mean quantum 
field. Such an interpretation is in concordance with the initial motivation of our analysis as 



directed to construct a covariant version of the chromo-magnetic models |T3 



Another conclusion which should be noticed is that the evaluated gluon mass determined 
from the allowed values of the condensate parameter g 2 C is tachyonic. The possibility for a 
tachyonic gluon mass has been recently discussed in the literature [JTSj, |26], |27J. This result 
indicates that the free propagator of the expansion is not reproducing itself after the first 
approximation. This property, however, does not mean a limitation of the approach. This 
is because, at variance with p2| |23|j , we are not assuming that the exact propagator has the 
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same delta function structure. This form only arose for the free propagator as determined 
by the Wick expansion around the considered initial state. 

The presentation is organized as follows: Section II is devoted to make precise the ele- 
ments characterizing the modified perturbation expansion. In the Section III the calculation 
of the (G 2 ) up to g 2 order for gluodynamics is presented and the result used to estimate 
the condensation parameter g 2 C. Further in Section IV the effects of the condensate on the 
quark masses are evaluated and the obtained values used to estimate the spectrum of the 



ground states in each baryon group as classified in [21]. In Section V the calculation of the 
effect of the condensate on the one loop self-energy for the gluons is evaluated. Finally, the 
results are summarized and some open questions and a conjecture are underlined in the last 
Section VI. 

II. THE MODIFIED FEYNMAN RULES 

In this introductory section, the main elements and conventions of the modified expansion 
introduced in [TB|, ^ [TJ| will be reviewed. The generating functional of the Green functions 
in the modified theory has the form 



Z[J,fj iVi lt] = i J 'D(^<r,c^tf r )exp{» J d'xC 



sources 



(1) 



C sources = £ + + -ra^a + -aja + ^rig + ^rt ^) 



N = J D (A r , c r , c r , f r , * r ) exp ji J d A xC^ 



(3) 



where sources for all the renormalized fields (A™, <T a , c ra , ty n , *ff n ) have been introduced 
in the usual manner, and the effective Lagrangian for the bare fields was selected in the 
Kugo-Ojima quantization procedure for gauge theory as 

£ = C G + C Gh + C Q (4) 

C G = ~F; V (x) (x) - -L (d»Al (x)) 2 , (5) 

C Gh = -td»c a (x)D; b (x)c b (x), (6) 

C Q = ¥ (x) {zYDv - m* Q 5«) ¥ (x) . (7) 

The sum over the six quark flavors will be omitted everywhere in order to simplify the 
exposition. We think that no confusion should arise for it. The gluon field intensity has the 
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usual form 



F; v (x) = d,A a u (x) - d u A« (x) + <f / a6c 4 (*) ^ (*) , 

where Df, Dj> are the covariant derivatives in the adjoint and fundamental representations 
of the group SU (3) respectively 

Df{x) = d^ ab -g*f ahc A^(x), 
D*j (x) = d^5 ij - tg*T ij ' a A^ 

Here, the bare fields and bare coupling constant and masses (<?*, Wq), were expressed in 
terms of their renormalized counterparts following the same conventions of Ref . |28| . Below, 
the renormalized masses and coupling will be denoted as tuq and g respectively. 

It should be recalled that as argued in |14], [H| the physical predictions for the value a = 1 
of the gauge parameter should have physical meaning whenever the adiabatic connection of 
the interaction does not lead the evolving state out the physical space. Therefore here we 
will fix this value of the parameter. 

Then, after the standard procedure of extracting out of the functional integral the ex- 
ponential of the terms of higher than second order in the fields (vertices) the generating 
functional takes the form to be employed in the evaluations below [[29|j : 



Z[j,fj,r),l£\ =iV- 1 exp| 



cG x3 qG x4 

D abc _j_ D abcd 



3H 3 5J a 5J b 5J c AH 4 5J a 5J b SJ c SJ d 

qGh r3 qQ 

U '-'in. 



Z [J,fj,7i,Z,Z], (8) 



with the free functional given by 

Zo [J, V, T], £, £] = exp ^i JaG G Jb + ifj r Gc h ris + 

and normalizes Z to be one at vanishing sources. Use will be made of the DeWitt compact 
notation p0 |, in which a Latin letter let say a,b,..., for a field symbolizes the space-time 
coordinates of it as well as all its internal quantum numbers. The same index in a source or 
a tensor symbolizes the same set of variables of the kind of fields associated with this specific 
index. For example A a = A a ^ a (x a ), c a = c a (x a ) and ty l = ^(xi). Note that a particular 
convention has been also employed in which the same letter is employed for the global index 
(for all the coordinated and internal indices) and the one indicating the internal coordinates 



in the explicit form of the quantities. Such a procedure was useful for the calculations 
and we also think that will not create confusion. As usual, repeated indices represent the 
corresponding space time integrals and contracted, Lorentz, spinor or color components. 
The following definitions also have been used in (|8[): 

5 5 5 



^ijk. 



5§i 5<£>j 5$ fc 



...S a ($) 



for a = G, Gh, Q and $ = A, c, c, \1> 
for a = G, Gh, Q, 



where G, Gh and Q mean the gluon, ghost and quark parts of the action respectively. 

There is only one main element determining the difference of the usual Perturbative QCD 
expansion and the modified one considered in JTB|, [TJ], |T5 |. It is related with the form of 
the gluon free propagator. As proposed in |T3[ there is an additional term to the gluon free 



propagator which is absent in the standard expansion. In [[[J, [L5| , such a term was shown to 
be a consequence of a Wick expansion based in a state constructed by acting on the usual 
vacuum with an exponential of pairs of zero momentum gluon and ghost creation operators. 
The ghost propagator |14], |15| could remain unmodified if the parameter left after fixing the 
form of the wave-function is taken as a real and positive one |14|, [15j . Such a selection will be 
employed here. The quark propagator is not affected in any way by the gluon condensation 
as introduced in |13|, [14|, However, as it will be commented in the summary, the results 
of the present work led to the idea about a possible physical relevance of introducing quark 
condensates along the similar lines as it was done for the gluons. We consider the exploration 
of this possibility to be exposed as one of the most interesting extensions of the work. 

A technical point can be noticed for precision. The operator quantization employed in 
HI, ffH] takes the ghost fields as satisfying the conjugation properties proper of the approach 



of Kugo and Ojima | 20| . However, at the level of the Feynman diagram expansion the 
difference with the standard procedure is only a change of variables. 



In accordance with above remarks and the results of |T3|, [14], |15[ the diagram technique 
defining the modified expansion has the following basic propagators for the renormalized 
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fields 



G ab = $ ab Q ,J "*^ b 
G j 



k 2 + is 



iC5 (k) 



® (m 2 — p 2 — is) 

and the same normal and counterterm vertices of the usual theory p8| . The dimensional 
regularization and minimal substraction (MS) is the renormalization procedure employed 
here. 

Being defined the Feynman rules and the conventions to be used, below they will be 
applied to the calculation of various quantities of interest. 

III. THE MEAN VALUE OF G 2 

The calculation of the mean value of the gluon field intensity squared is considered in 
this section. It extends the former evaluation done in |T3[ within a simpler approximation. 



The result will be employed afterwards for fixing the value of the condensation parameter 
g 2 C to be the necessary one for furnishing the currently estimated value for (g 2 G 2 ). The 
calculation will be performed up to order g 2 including all the dependence on the condensate. 

It is a helpful circumstance here, that the one loop renormalization procedure of the 
modified expansion can be performed identically as in the standard PQCD. This is a direct 
consequence of the fact that the (^-functions kill the only existing integrals at the one loop 
level and therefore all the corresponding counterterms in the dimensional regularization are 
identical. This fact will be helpful in regularizing and renormalizing the composite operator 
(G 2 ). 

For the evaluation of (G 2 ) the following definition of this quantity will be employed 

2 = (0|g g |0> = / D (A r , tf, g, Or, vF) S g [AT\ exp {z / d'xC*— } 
[ ' fdx D Jdx D f D(A r ,c r iC r ,V r ^ r )exp{iJd 4 x£} ' U 

where S g [A r ] represents the gluon part of the action in the absence of the gauge breaking 
term depending on the gauge parameter a, evaluated at the renormalized gluon fields A r , 
that is 



S g [A r \ = -\ J d*xF™(x)F r ^ a (x) 



Therefore, (G 2 ) will be considered as defined by the mean value of the gluon action 
divided by the space-time volume in dimensional regularization. 

After substituting in ([5]) the perturbative expression for the mean value in terms of the 
Generating Functional (H) it follows 



<0| S 9 |0> 



C9 X2 Q9 X3 q9 
°afe , °afec |_ °abcd 



2i 2 8J a 8J b 3\i 3 5JJJ b 8J c 4H 4 6JJJ b 5J c 5J d 



z [J, v, v, 6, 6] 



(10) 



In the momentum representation the kernel S 9 ab has the form 

S° b = (2vr) 4 5 ft, + k b ) 5 ab [- (g^k* - K,k a ^)] 



and the 5f bc and S 9 abcd ones are the usual three and four legs gluon vertices . 

Let us consider below the evaluation of (G 2 ) up to order g 2 . As any Feynman diagram 
can be expressed as a polynomial in the constant C, after associating a new diagrammatic 
line for this part of the free propagator, all the Feynman diagrams can be divided in two 
groups. The first, let say G\ including the diagrams not showing the new propagator line 
and another group, G2 having one or more, let us call them, "condensate" lines. Such lines 
have associated to them a multiplicative Dirac delta function at zero momentum. 

After writing out all the expressions of type G2 (up to second order in the renormalized 
coupling g 2 ) contributing to (0); it follows that the only non- vanishing diagram is the one 
appearing in the Figure [I] in which both gluon lines are substituted by "condensate" lines. 
The remaining terms vanish thanks to the zero momentum evaluation produced by the 
Dirac Delta functions introduced by the "condensate" lines and also due to the vanishing 
of tadpoles in dimensional regularization. Also, it should be underlined that the retention 
of the ie factors defining the Feynman propagator, up to the evaluation of the integrals was 
also important in obtaining the zero result for some of the diagrams in G2- 




FIG. 1: Diagram contributing to G2 terms in the gluon condensate parameter. 
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The contribution related with Figure p] had been calculated before and the result is |L3|, [14 

_128% 2 C 2 

Next, the diagrams in G\ are the ones contributing to G 2 up to order g 2 in standard 
perturbative QCD. They are illustrated in Figure |2|. The terms T$ and T-j represent the only 
non-vanishing and non-tadpole like second order contributions. It can be noticed that the 
appearance of two tadpoles in the term T-j leads to the exact vanishing of this contribution 
in dimensional regularization. 





•OO 



6 A 7 

FIG. 2: Diagram contributing to G\ terms in the gluon condensate parameter. 



For the diagram T$ corresponds the expression 



j_ r d P D <r<?* (fl) 

Te 3! i 3 J (2tt) d z p 2 + is ^ 



ab 



(P), 



where is the gluonic contribution to the polarization tensor in second order, which has 
been explicitly calculated. After substituting the known result for II and performing the 
Wick rotation, the previous expression turns out to be of the form 

dq D 



E' 



(27T) 



D 



(Q 2 



(11) 



Then considering the formula [31| 
dq D 1 1 i{-l)- a - l3 {-k 2 y a ~^ 2 T{a + (3 -n/2) 



(2tt) d (q 2 ) a (q - k) 2f3 
B(n, m) 



(16vr 2 ) n / 4 



D{--(3,--a). 



r(m)r(n) 
T(m + n) 
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for the case a + (3 = 2 — y and fixing — 0, it follows the vanishing of diagram T§ within 
the dimensional regularization for D > 2. 

It rest now to discuss the contributions of the tadpole like diagrams. It is worth noticing 
at this point that they are quite similar in their form to the one corresponding to the 
square of the electromagnetic field strength in QED. In fact they differ from it only in some 
numerical color factors. An important quantitative difference between the expressions for 
both theories is related with magnitude of the coupling constant which is very much higher 
in QCD. Let us consider the T\ — T 5 terms in Figure [2] embodying the gluon, ghost and 
fermion contributions with the inclusion of the one loop counterterms rendering finite the 
polarization operator at one loop This expression turns to be of the form 



d P D (p^Pu-p 2 g^)(g^ 



pfp" 



J (ZTTj^ p A + is 

where II as usual is defined through the polarization tensor 

Ilf u (p)=6 ab (p fl p iy -p 2 g, 1/ )U(p 2 ), 

and it can be expressed as the superposition of the pure gluodynamical contribution and 
the quark one as |28[ 



n( P 2 )=n G ( J 9 2 ) + n Q ( F 




-(2-D/2)n 



+(^3-1), (13) 

where Cg = 3 and Tr = 1/2 for SU(3) and /i is the scale mass in the dimensional regular- 
ization. 

Considering first the pure gluonic contributions associated (1 — II G (p 2 )) in (|l|)> it follows 
that they reduce to a sum of dimensional regularized integrals of the same form as ((T/TJ) . Thus, 
the second order gluodynamical contribution of the tadpole like contributions vanish. Let us 
analyze now the quark contributions associated to LTq. It can be noticed that for massless 
quarks, that is rriQ = 0, the LTq has the same analytical dependence of p 2 as the gluonic part 
Ug and all the quark contributions vanish after the regularization is employed. However, for 
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massive quarks the result for the second order dimensionally regularized contribution along 
the same lines employed for the previous terms gives rise to a non-vanishing contribution 
of the following form after the substraction of the divergent pole terms at D = 4. The 
expression for the regularized term before substraction is 



9V' d e T6 = £Kr 2 



6A^g^ 8 - 2D T R (D - l)r(2 - D)T(4 - D) 



2 7 ~ 2 - D (47r) £ T(! — D) 

This expression has simple as well as double poles which after substracted (in the D 



4 limit following the rules of the MS scheme for composite operators [32]) give a finite 
correction to (g 2 G 2 ) for the gluon condensate parameter. It can be evaluated and the 
numerical expression becomes 



(gY- D E T6 ) fin = 0.000642V ( J>Q 



55.32 - 37. J In | A- | +6.0 In 2 [ 

m% \ mi. 



However, it should be noticed that this finite expression has been obtained for the sec- 
ond order correction of the polarization operator. The summation of the infinite ladder 
of self-energy insertions leading to the one loop propagator can be suspected to produce a 
noticeable modification of this result. After also considering the high indefiniteness in the 
quark current masses, and the technical difficulty we had in evaluating the dimensionally 
regularized expression for the above mentioned self-energy insertions, we decided to employ 
for the further consideration here the purely gluodynamical contribution. In support of this 
way of proceeding there is also another circumstance. It is needed to notice that a refined 
version of the definition of the composite operator G 2 (x) could be suspected to exist in 
such a way that the terms associated to a tadpole diagram of the exact propagator vanish. 
This property would assure the exact vanishing of the (F 2 ) in QED a property that to our 
knowledge, is expected. Because, at one loop, both QED and QCD tadpole expressions 
are basically the same, the above written quark contribution would not really appear. The 
existence of the mentioned definition, opens the possibility for the second order gluodynam- 
ical contribution to be exact one for QCD in this approximation. However, we shift the 
discussion of this question to further analysis to be considered elsewhere and simply adopt 
the gluodynamical result. 

Therefore, the total second contribution to the mean value of G 2 in gluodynamics is 
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defined by the only non- vanishing term 

2 _ 288g 2 C 2 

coming from the condensate contributions arising from the diagram Tq which is defined by 
the four gluon interactions. Finally, let us determine in a phenomenological way a value for 
the constant g 2 C by making use of a currently estimated value for (g 2 G 2 ): 

(g 2 G 2 )^ 0.5 (GeV/c 2 )\ 

Therefore, the parameter g 2 C takes the value: 

^ 2 C = 64.9394 (GeV/c 2 ) 2 . (14) 

This fixing of the condensation parameter allows to investigate the physical predictions 
of the modified expansion in the next section. 



IV. CONSTITUENT QUARKS MASSES 

Let us consider now the effect on the current quark masses produced by the condensate. 
The full inverse propagator for the quarks in terms of the corresponding mass operator has 
the general expression 

G Q 1 (p) = (m Q -p^-J:(p)), (15) 

where the color index identity matrix has been omitted and £ (p) is the mass operator. Its 
expression is related with the diagram in Figure |3|. 

i} u 

FIG. 3: Diagram for S (p) in one loop. 

The form for E (p) up to order g 2 , within the modified Feynman rules can be written as 
E (P) = 9 2 Cf f d \ ^^Q + ^-^^^g W, (16) 

in which the color factor Cf for the SU (N) group is Cf = which for the case of 

interest here SU (3), reduces to Cf = I [p8fl . 
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Taking into account the modified gluon propagator and the standard relations for the 
7 matrices (which here are adopted in the same way as in Ref. pBf ) the one loop mass 
operator expression is simplified to be 

J (2vr) i (m 2 Q -ip-ky) \k 2 J 

Let us disregard for the moment the term not involving the condensate parameter C in order 
to first study the sole influence on the quark masses of the condensation effect. Then, the 
expression for E reduces to 

g 2 C F C2{2m Q -p a la ) 
[) (2vr) 4 (m 2 Q -p 2 ) ' 
Therefore, the inverse propagator for the quark takes the simple form 

( M 2 \ ( M 2 \ 

where the constant M is defined by 

m2 = 2g 2 C f C = UGeV/c 2 ) 2 , 
(2tt) 4 V ' ' 



in which the numerical value has been obtained by using (|14|). 

The zeros of the determinant associated to the inverse propagator (|I7D then allows to 
determine the effects of the condensate (reflected by the parameter M), on the effective mass 
of the quarks. Let below s,s — 1...6 an index indicating each kind of quark characterized 
by its particular current mass rriQ s . 

In order to find the dependence of the effective quark mass m 2 . as a function of the 
current mass parameters vhq s an analysis was done for the zeros of the determinant of the 
inverse green function (0). It follows that there is only one solution having a squared mass 
being positive for arbitrary values of the current quark masses. The existence of this real 
and positive solution for the squared mass becomes possible thanks to the real and positive 
character of C, as was shown previously JL4|, jl5| . There is another interesting branch of 
the zeros of the determinant (dispersion relations) showing vanishing values for the effective 
mass in the light current mass limit. For the purely real solution the value of the effective 
quark mass m qs = f (rriQ s ) as a function of mq a is shown in the Figure ^. The graph is 
plotted for the region rriQ s < 2 GeV/ c 2 which contains the current mass values of the u, d, s 
and c quarks. 
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FIG. 4: Real solution for the quark mass as a function of the Lagrangian mass (masses in GeV/c 2 ). 



As it can be appreciated in the picture, the calculated effective quark masses for light 
flavors (u, d and s) are clearly predicting the values of the quark masses being in use in the 
constituent quark models of the baryons. That is, the light quarks get a weight of near one 
third of the nucleon mass. In the table [I] the mass values obtained for each of the quark 
flavors are shown pil. From the global properties of table | it can be observed that the main 
effect of the gluon condensate seems to be the dressing the light quarks with a cloud gluons 
having a total mass of one third of the nucleon mass. Then, the results point in the direction 
of the mainly glue nature of the constituent quark masses of the u, d and s quarks within 
many baryon resonance. These results then support the idea that a modified perturbative 
expansion like the one being considered, in which the effect of the gluon condensate have 
been incorporated, would be able to predict with reasonably good approximation properties 
of the low energy strong interactions. Similar values for the constituent quarks masses have 



been obtained by different methods in |33|, |34 



We will now investigate the influence on the mass spectrum produced by including the 



standard one loop contribution in the dispersion relations. After substituting the known [31 



quark one-loop self-energy expression, the one loop inverse quark propagator can be written 
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TABLE I: Quark mass values in presence of the condensate in units of MeV/c 2 . 



Quarks 


Exp a 
m Low 


Exp b 
m Up 


^Theo c 
m Med 


(u) 


1.5 


5 


334.944 


(d) 


3 


9 


336.287 


(s) 


60 


170 


388.191 


(c) 


1100 


1400 


1315.241 


(b) 


4100 


4400 


4269.572 


(t) 


168600 


179000 


173800.48 



"Reported lower bound value for the Lagrangian mass 
b Reported upper bound value for the Lagrangian mass 

c Calculated mean value of the constituent mass, for the lower and upper bound of the Lagrangian masses 
as determined only by the condensate contribution to the self-energy in the one loop approximation 



as 



Gq 1 (p) = (m Q - - E (p)) = A{p 2 ) m Q - B{p 2 ) p" 7/i 



(18) 



m Q 



1 + 



AC F g 2 C 



(2tt) 4 m 2 Q -p 2 



{An) 2 



x 



-6 + 47 - ln(4vr) + 4 In 



2 2 

m Q~P 



-4 



m 



Q 



In 



P 



2 2 
m -V 



1717 



1 + 



2C F g 2 C 



(2tt) 4 m 2 -p 



1 - 7 + ln(47r) + - In 



(4tt) 2 
p 2 ' 



+ 



o o / m n — P 
-) In ' Q ^ 



p^ \ /i" / p^ 

where the dimensionless coupling constant g 2 has been substituted by its running expression 
as a function of the scale parameter \i as 

2 _ 4vr 0.12 

9 ^> ~ i + To^j r-jj-y- 

x ~ 4tt 111 V 91.18/ 

That is, the value of a 2 (/i) = (? 2 (/i)/(47r) is fixed to be 0.12 at the scale of the Z mass 
/i = m z = 91.18 GeV. In order to automatically fix the scale /x to the one associated to the 
magnitude of the obtained solution for the mass, the value of /1 was taken to satisfy /1 2 = p 2 . 
That is, to reproduce the same value of the solution for the squared mass. 
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Before proceeding, it should be noticed that the inverse propagator, after the introduction 
of the one loop terms, gets a branch cut at the real axes for the considered complex variable 
p 2 for Refjo 2 ) > uiq. However, the ie prescription selects one of the sides of the cut as giving 
a well defined value of the inverse propagator at real p 2 values. 

After multiplying ( |ISD by A(p 2 ) mq — B(p 2 ) p^ r y ll the modified mass shells can be most 
easily found by solving 

D(q 2 ) = A 2 (p 2 ) m 2 Q - B 2 (p 2 ) p 2 = 0. 

Concretely, we estimated the corrected values of the masses in Table [I] by searching for 
the peaks of the inverse for the absolute value of D(q 2 ) as a 2D-function of the real and 
imaginary parts of p 2 . The new values are reported in Table O. As it can be observed the 
masses were not drastically changed. They all show a relatively small reduction in their 
magnitude, then indicating the consistency of the predicted constituent masses for the light 
quarks at the one loop level. 

TABLE II: Quark mass values in presence of the condensate in units of MeV/c 2 . 



Quarks 


Exp a 
m Low 


Exp b 
m Up 


Theo c 
" L Med 


(u) 


1.5 


5 


309.15 


(d) 


3 


9 


309.15 


(s) 


60 


170 


354.00 


(c) 


1100 


1400 


1343-200f 


(b) 


4100 


4400 


4515-65i 


(t) 


168600 


179000 


190000-1450i 



"Reported lower bound value for the Lagrangian mass 
^Reported upper bound value for the Lagrangian mass 

c Calculated mean value of the constituent mass, for the lower and upper bound of the Lagrangian masses 
after the inclusion of the standard one loop correction for the self-energy. 

An interesting outcome following from the inclusion of the standard one loop terms, 
should be stressed. It consists, in that for the smallest rriQ masses associated to the lower 
and higher bounds for the u and d quarks, \ D ( q z^ although showing a definite peak at 
0.30915 GeV does not really diverge near this value. Therefore, the propagator is not 
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showing a definite pole associated to a propagating quark mode at this p 2 value. However, 
at the level of the Bethe-Salpeter equation for bound states of quarks, it should be expected 
that p 2 at this peak should behave as the non-relativistic value for the square of the mass 

2 

entering the kinetic energy expression Henceforth, in spite of the presence of free 

propagating modes in the simplest approximation, the inclusion of the quantum one loop 
correction makes the u and d modes non-propagating (confined) in the same way as it occurs 
in the approach of Munczek and Nemirovsky |22|, . 

For greater values of mq in the range presently estimated for the s quark current mass 
(0.060 — 0.170 GeV) in the neighborhood of the non-diverging peak of r^^ji , this quan- 
tity starts showing divergency points at complex values of p 2 which represent damped but 
propagating quark modes. The just described structure of the spectrum suggests a possible 
explanation for the drastic difference of stability between the u and d composed hadrons 
and the resonances including s, c, b and t quarks. The point is that the presence of poles 
at complex p 2 for the heavier quarks could lead to the existence of disintegration channels 
for the higher mass resonances in which the s, c, b and t quarks can escape in a damped 
propagating mode to be after hadronized by some separate mechanism. On another hand, 
as the u and d quarks do not show propagating states, that could remain in a perennial 
confinement within an undisturbed nucleon. A further study of the consistency of this ex- 
planation is expected to be done within the planned extension of this work to investigate 
bound state equations for hadrons. 

It can be concluded that the constituent quark masses evaluated from the simplest ap- 
proximation for the dispersion relation considered above are stable under the inclusion of 
the standard one loop corrections to the self-energy. 

Next, in the Table |H| the masses of the ground state resonance selected each of them 



from one of the groups of baryons as classified in [pi]] , are shown in comparison with the 
values obtained for them by a theoretical estimate consisting in adding the here calculated 
masses for each of their known quark constituents. The quark masses of the Table | were 



employed for this evaluation. Also, in the Table IV, the masses for some vector mesons 



reported also in are shown in conjunction with the result of the sum of the calculated 
masses for their constituent quarks (Table |). 

As it can be appreciated the results obtained for the vector mesons and baryon resonance 
are only estimates. However, the obtained values are in a reasonable well correspondence 



19 



TABLE III: Experimental and Theoretical Baryonic Resonance Masses in units of MeV/c 2 . 



Baryon 


Exp.Val. 


Th.Mean.Val. 


Rel.Err. 


p(uud) 


938.27231 


1006.175 


7.24 


n(udd) 


939.56563 


1007.519 


7.23 


A(uds) 


1115.683 


1059.422 


5.04 


£+(uus) 


1189.37 


1058.078 


11.04 


£°(uds) 


1192.642 


1059.422 


11.17 


E-(dds) 


1197.449 


1060.766 


11.41 


S°(uss) 


1314.9 


1111.325 


15.48 


E-(dss) 


1321.32 


1112.669 


15.79 


fl~ (sss) 


1642.45 


1164.572 


29.10 


A+(udc) 


2284.9 


1986.472 


13.07 


H+(usc) 


2465.6 


2038.375 


17.33 


^(dsc) 


2470.3 


2039.719 


17.43 


^O(ssc) 


2704 


2091.622 


22.65 


A^(udb) 


5624 


4940.803 


12.15 



TABLE IV: Experimental and Theoretical Masses for a group of Vector Mesons in units of MeV/ c 2 . 



Meson 


Exp.Val. 


Th.Mean.Val. 


Rel.Err. 


„ / iiu—dd \ 
P \ V2 ) 


770.0 


671.231 


12.83 




781.94 


671.231 


14.16 


(f>(ss) 


1019.413 


776.381 


23.84 


J/ii>(lS) (cc) 


3096.88 


2630.482 


15.06 


] Y(lS){bb) 


9460.37 


8539.144 


9.74 



with the reported experimental values. It can be also noticed that the validity of the present 
calculations implies that the binding energy contribution to the baryon rest masses should 
be small. This property then could justify a non relativistic description of such resonances. 
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As it was remarked before, another form of the dependence of the effective quark mass 
with rriQ s is also following from the dispersion relations. It, by contrast shows a vanishing 
effective mass when mq s — > 0. The existence of this branch opens the possibility for the 
existence an appropriate description of the light mesons (II and K) within the proposed 
scheme. They could be associated to bound states of light quarks propagating in such 
modes. The analysis of this question deserves a separate treatment to be examined elsewhere. 
The graphic for the mentioned branch is represented in Figure for the region of rriQ s < 
0.1 GeV/c 2 , a point at which the change in the analytic dependence signals where the 
solution gets an imaginary part when the current mass value itiq grows. 



0.25 



0.2 



0.05 




0.02 0.04 0.06 0.08 0.1 

mQ 



FIG. 5: A solution for the quark mass as a function of the Lagrangian mass, that satisfies m qs — > 
when rriQ s — > (masses in units of GeV/c 2 ). 



V. TACHYONIC GLUON SELF-ENERGY 

Let us consider now the evaluation of the gluon self-energy. As it is known the exact 
gluon propagator can be written in the form 

sab / n - p ^ p n \ 

GKyj p 2 I i + n ( P 2 ) p 2 J 
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where the function IT is defined at the one loop level, by the expression 



which is diagrammatically represented in Figure ^. 




a ,' 




FIG. 6: Diagrams for U ab (p 2 ). 

Again, expanding in powers of the condensate parameter C, the total contribution can 
be decomposed in the standard one (associated to the ghost and the quark loop integrals) 
and the C dependent ones Hq and IT^ which can be readily calculated to be 



jrab _ 9 2 CN , / p P/iaPwA 



(2tt 
(2tt) 

Summing up these particular two terms leads up to 

^{v 2 )=^{p, a p, b -g^ b p 2 )^- 2 . 

(2%) p 2 

Considering again our case the SU (3) group, that is N = 3, and disregarding in a first 
consideration the standard one loop contribution, the function II reduces to 

6g 2 C 1 

which in turns leads out to the following condition for the poles of ( |19|) 

p 2 -m 2 G = 0, 

where 

Q9 2 C _ nnK ,^„*s2 



n (p 2 ) = ( 2 o) 



m 



G ~ (2vr) Z 



-0.25 (GeV/c 2 )' 
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Therefore, as the parameter C has been defined before as real and positive [14], [15[ , 
it follows that the transverse gluon mass correction becomes tachyonic. The ability of a 
tachyonic mass in producing improvements in models for the inter quark potential has been 
recently argued in the literature B7[|. The effect is linked with the introduction of a 
linearly rising term as a first correction to the Coulomb potential in the massless case J27 



Therefore, the tachyonic result arising here for the gluon mass appears to be of interest for 
a derivation, within the present framework, of former successful phenomenological bound 
state models for mesons |35|, |37H . We are grateful to J. P. Vary and J. Qiu for the 
suggestion of this possibility. It should be also underlined that the tachyonic character 
of the gluon self-energy arising here, is fully compatible with the more general result of 



Fukuda for the normal Green functions in QCD, when gluon condensation is present |38 
This circumstance then suggests, that the modified expansion under discussion could be a 
sort of perturbative representation of the Fukuda non-perturbative treatment of the gluon 
condensation in QCD. In a different manner, the approach of Celenza and Shakin, and 
other authors can be expected to be related with a different initial state also incorporating 
gluon condensation. It could be the point that, at the end, both pictures can appear to be 
equivalent at the level of exact results. This could be the case because the physical degrees 
of freedom are not the gluons and quarks ones, but their colorless composite combinations 
which physics is determined by the bound state dispersion relations. 

Finally in this section it deserves to be stressed that if a modification of the ghost prop- 
agator is present, then a mass terms for the longitudinal propagator could also appear and 
the polarization operator (gluon self-energy) becomes non-transverse. This result means a 
violation of a Ward identity representing the gauge invariance by the presence of ghost (and 
fermion) condensates. Thus, it follows that the selection of the parameter of the initial state 
in |L4], [15|] was the appropriate one for the satisfaction of a gauge invariance condition at 
this stage. The physical meaning of the alternative procedures of constructing the initial 
state needs however for clarification. 

VI. SUMMARY 

The following results can be enumerated: 

1) A calculation of the mean value of G 2 up to order g 2 in gluodynamics, was performed. 
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The result coincided with a previous one obtained in a simpler approximation [13, [14| and 
which signals that quartic gluon self-interaction is relevant in fixing a non-vanishing gluon 
condensation parameter already in the first order of the perturbation theory |[19| . 

3) The one loop corrections for the quark masses under the influence of the condensate 
were evaluated. When only the terms reflecting the effect of the condensate are considered, 
the results for the constituent masses showed a branch which is real as a function of all the 
values of the quark current mass. Another solution which is real up to a critical value of 
the Lagrangian masses, and having a vanishing proper mass in the limit of small current 
masses is also obtained. The former branch predicts values of one third of the nucleon mass 
for the masses of the light quarks. The effective quark masses for all the quark flavors were 
evaluated as determined by their known current mass in the mentioned approximation. 

4) After the inclusion of the rest of the terms completing the full one loop contribution to 
the quark self-energy, an interesting outcome arises. The u and d quarks, although showing 
high peak values of the propagator components near the former values of the squared proper 
masses, become confined. That is, their propagators do not show pole structures. On the 
other hand the s, c, b and t quarks acquire poles but now not at real values of p 2 . Therefore, 
these free modes are damped. These properties suggest a possibility for the explanation of 
the disparity of stability properties in the nucleons and higher resonances. It would work 
as follows, the absence of propagation solutions for the u and d quarks could assure the 
unlimited life time of nucleons. Moreover, the damped oscillations present for the rest of 
the quarks can create decay modes for the bound states composed of them. 

5) By mean of the obtained mass values, the ground state energies selected within each of 



the groups hadron resonance reported in |2T| , were estimated through the simple addition of 
the masses of their known quark constituents. It could be remarked that the obtained values 
are reasonable ones in comparison with the data, after taking into account the high errors 
in the experimentally determined Lagrangian masses. This fact supports the hypothesis 
about the main glue nature of the constituent mass values of the light quarks. The small 
contribution of the binding potential energy of such resonance to their total mass is also 
indicated and thus suggests the applicability of non relativistic approximations in their study. 

6) The contribution of the condensate to the gluon self-energy of order g 2 was calculated. 
A tachyonic result for the gluon mass arises for the parameter values employed in the con- 



struction of the initial state defined in p4], |T5||. This conclusion is in concordance with the 
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analysis of Fukuda for gluon condensation f38| . As it was remarked before, there are recent 



arguments claiming the possibility of a tachyonic gluon mass |26| , p7[ and its role in 
improving the quark antiquark potential in bound state models of mesons p? . 



It seems useful to comment here some basic issues related with the modified expansion 
being considered: 

a) The Feynman expansion depending on the condensate parameter and taken for the 
gauge parameter a = 1 corresponds with the Wick expansion around a physical state of the 
interaction free QCD. In this somewhat limited way the arbitrariness of the gauge parameter 
in the result for the gluon mass in [|l^| was solved in (14], |15| . However, in order to furnish 
a fully consistent picture including renormalization, a physically equivalent version of the 
expansion for arbitrary values of a should developed. This question will be addressed in 
future works. 



b) In fi"4| , |15|| the condensation parameter C was defined as a real and positive one. This 
result determined the tachyonic character for the gluon mass found here. Further, the sign 
of C also fixes the real values of the order of the one third of the nucleon mass for the 
constituent quark masses also calculated in the present work. The opposite sign produces 
confined quark modes. 

c) An aspect which is important to underline is that the selection of the parameters 
defining the initial state in |TJ|, ^] was designed also in order to impose the absence of 
a modification for the free ghost propagator. However, as it followed from the present 
analysis, this property in turns, is related with the fact that a condensate modification of 
the ghost propagator can produce a longitudinal contribution to the self-energy. But, such 
a term should not exist for the transversality condition of the polarization operator (a Ward 
identity) to be obeyed. Thus, its appearance could break manifestly the gauge invariance. 
Therefore, the present work is also given foundation to the non-modified ghost propagator 



choice considered in |L3], [14| p~5| . 

At the present stage we identify some important questions to be addressed in the further 
extension of the work: 

1) To explore the possibility of showing the exact gauge parameter independence of the 
physical results within the modified expansion. 

2) To address the general proof of the renormalizability in the expansion. 

3) To search for a derivation of existing successful bounded state models for the heavy 
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quark mesons |36|, |37l , as realized by the a ladder approximation for the Bethe-Salpeter 
equation within the proposed modified expansion. The presence of a tachyonic gluon prop- 
agator in the approach (which is argued to have the effect of introducing a linearly rising 
component in the inter-quark potential |2"7|| ) and the obtained constituent values for the 
light quark masses already support such a possibility. 

4) Improve the study of the effective action as a functional of the condensate parameter 
done in ]13| by also including the gauge field in order to search for a variant of the leading 
logarithm model useful for the investigation of field configurations associated to inter-quark 
strings, nucleon-nucleon potentials, etc. 

5) It can be underlined that the considered expansion seems that could furnish a general 
way for the perturbative description of pair condensates in Many Body theory. An important 
example coming to the mind is the BCS state of the usual superconductivity. Pair states 



being more closely connected with QCD were discussed in [^, |4(|. In addition, in [g(J it 
was advanced before the idea of using diquark condensates as generating an alternative for 
the Higgs mechanism. 

Finally, we would like to conjecture about a possibility suggested by the results of this 
work. It is related with the question about whether the mass spectrum of the whole three 
generations of fundamental fermions could be predicted by a slight generalization of the 
modified perturbation expansion under consideration. In this sense, the given arguments 
led us to the idea that (after the introduction of quark condensates along the same lines 



as it was done for gluon ones in |13[ [14], |i5| ) the received perturbation expansion can have 
the chance to predict both, the Lagrangian mass and the constituent quark mass spectra 
of the three families of fundamental fermions. The fermion condensates as described in 
the proposed perturbative way, would have the role of producing the Lagrangian quark 
masses, through the chiral symmetry breaking. The gluonic condensates, in one hand, and 
as illustrated here, could be responsible of generating states of large constituent mass for the 
low mass quarks states. In another hand, it seems feasible that the higher order radiative 
corrections (including color interactions with the condensate as well as chiral symmetry 
corrections) could also determine the mass spectra for leptons and neutrinos. The smaller 
scale for the masses of these particles could be produced by the lack of lower order either 
color or electric-weak interaction terms in their self-energy. Therefore, the possibility that 
the Lagrangian mass spectrum of the three generations of the fundamental fermions could 
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be predicted by a modified perturbation expansion of the sort being proposed is suggested. 
Work directed to investigate this possibility will be considered elsewhere. The modified 
expansion under study could be a useful technique of implementing this idea. 
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